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3 Using particle-in-cell simulations and kinetic theory, we explore the current- 

4 driven turbulence and associated electron heating in the dissipation region 

5 during 3D magnetic reconnection with a guide field. At late time the tur- 

6 bulence is dominated by the Buneman and lower hybrid instabilities. Both 

7 produce electron holes that co-exist but have very different propagation speeds, 
a The associated scattering of electrons by the holes enhances electron heat- 

9 ing in the dissipation region. 
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1. Introduction 

10 Magnetic reconnection is the driver of explosive events in nature, such as solar flares, 

11 substorms in the magnetosphere of the Earth and flares from magnetars and the accretion 

12 disks of black holes. Satellite observations in the Earth's magnetosphere indicate that 

13 magnetic reconnection drives turbulence. Electron holes, which are localized, positive- 
« potential structures caused by plasma kinetic instabilities, have been linked to current 
is sheets associated with magnetic reconnection in the magnetotail [Farrell et al, 2002; 

16 Cattell et al, 2005; Andersson et al, 2009], the magnetopause [Mafcwmofo et al, 2003], 

17 and the laboratory [Fox et al, 2008]. Lower hybrid (LH) waves and other plasma waves 
is appear in conjunction with electron holes in the magnetotail events. Electron holes can 
w scatter electrons, causing heating and possibly anomalous resistivity to facilitate fast 

20 magnetic reconnection. 

21 During magnetic reconnection, a parallel electric field generated around the x-line drives 

22 electron beams. Simulations with a guide field show that these intense beams can drive the 

23 Buneman instability, which forms bipolar structures in the parallel electric field [Drake 

24 et al, 2003]. Later in time transverse electric fields develop. Following a suggestion 

25 that these transverse fields were current- driven lower hybrid waves (LHI) [McMillan and 

26 Cairns, 2006], Che et al. [Che et al, 2009] showed that both the LH and electron- 

27 electron two-stream instabilities resonate with the high velocity electrons and therefore 

28 dominate the interactions with the highest velocity electrons in narrow current layers. 

29 Which instabilities develop during reconnection and how they interact remains unknown. 
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30 During magnetic reconnection we demonstrate that two distinct classes of electron holes 

31 with very different propagation speeds exist simultaneously. Slow moving holes are driven 

32 by the Buneman instability and at the same time and locations fast moving holes are 

33 driven by the LHI. Both take the form of nonlinear Bernstein-Greene-Kruskal (BGK) 

34 solutions [Bernstein et al, 1957] since the measured bounce time of electrons in the holes 

35 is short compared with the hole lifetime. The trapping and scattering of electrons by 

36 holes of disparate phase speed enhances dissipation during reconnection. 

2. Simulation 

37 We carry out 3D magnetic reconnection simulations with a strong guide field similar to 

38 those carried out earlier [Drake et al, 2003] but with a much larger simulation domain: 

39 L x = 4di, L y = 2di, and L z = 4di, where di = c/u pi and ui p j is the plasma frequency of 
M a particle species j. The reconnecting magnetic field is B x /B = tanh[(y — L y /A)/wo] — 

41 tanh[(y— 3L y /A)/wo\ — 1, where B is the asymptotic amplitude of B x outside of the current 

42 layer, and w is the half-width of the initial current sheet. The guide field B 2 = B 2 — B 2 is 

43 chosen so that the total field B is constant. In our simulation, B is taken as 26 1 / 2 B . The 

44 initial temperature is T e = Tj = 0.04rriiC 2 A , the ion to electron mass ratio is 100, the speed 

45 of light c is 20ca with ca = Bo/(Annom, i ) 1 / 2 , the Alfven speed. The initial drift speed of 

46 4ca is just above the electron thermal speed 3ca and marginally exceeds the threshold to 

47 trigger the Buneman instability. 

48 Magnetic reconnection induces a parallel electric field around the x-line and drives 

49 an intense electron beam. At flit = 3 (fli = eBo/rriic), the electron beams have been 
so accelerated to IOca and to IAca at flit = 4. We show the current sheet around the x-line 
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51 in the x—y plane at Qrf = 3.3 in Fig. 1 (a). At the beginning of the magnetic reconnection 

52 simulation, the Buneman instability with wavevector along the magnetic field z direction 

53 is excited. In the cold plasma limit, the phase speed is (m e / (2m i )) 1 / 3 \v c i z \/2 ~ lc A and 

54 the growth rate is 7 ~ v / 3w pe (m e /(2mj)) 1 / 3 /2 ~ 29f2j [Galeev and Sagdeev, 1984]. The 

55 Buneman instability saturates within a short time. Later in time two distinct spatial 

56 structures of the electric field are observed: localized bipolar structures dominate E z 

57 and long oblique stripes dominate E x . A surprise is that there are two types of bipolar 

58 structures. At 0,$ = 3 one has a velocity close to zero and the other moves with a velocity 

59 of 3ca- By Hit = 4 the velocity of the second increases to 7ca- In Fig. 1 (b, c) we show 

60 E z and E x in the midplane x — z of the current sheet at VL(t = 3.3. The structures move 

61 to the left in this figure, which is in the direction of the electron drift. The downward 

62 (upward) arrows point to fast (slow) moving electron holes. To see the two classes of 

63 holes more clearly, in Fig. 2 (a, b) we stack cuts of E x (z) and E z (z) at the x-line versus 

64 time. The dark and light bands mark the development of the bipolar structures seen in 

65 Fig. 1 (b, c). The slopes of these bands are the phase speeds of the waves. During the 

66 time interval Qit = — 2, the phase speed of the waves increases, which was expected 

67 since the streaming velocity of the electrons increased as the reconnection driven current 

68 layer shown in Fig. 1 (a) developed. During the time interval fl^t = 2 — 4 two distinct 

69 phase speeds, particularly in E z , are evident. In Fig. 2 (b) the structures cross each 

70 other at the same value of z, which indicates that this result is not due to the spatial 

71 structure of the streaming velocity. In Fig. 3 we show E z and the z — v z phase space 

72 around (x,y) = (1.2dj, 1.5dj) at £1$ = 3 to reveal the structure of the fast moving holes. 
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73 There are no slow holes in this region at this time. In (a) the most intense hole is marked 

7 4 by the arrow. In (b) the center of the z — v ez phase space of this bipolar structure is 

75 marked by the star. The electrons encircling the star indicate that electrons are trapped 

76 by the bipolar field. The strong electron heating due trapping is evident. 

77 Electron holes in the simulation exhibit a complex dynmics: formation, dissipation and 

78 reformation. The lifetimes 77 of the two classes of electron holes are distinct, around O.lfi" 1 

79 and 0.2J1" 1 for the fast and slow holes, respectively. In both 77 exceeds the bounce time 
ao of the trapped electrons, r b ps y/m e X b /y/2e5E z ~ 0.02Q" 1 , where \ b is the characteristic 
si wavelength of the electron hole. Thus, electron trapping takes place and we therefore 

82 interpret the holes as BGK structures. [Bernstein et al, 1957]. 

3. Kinetic Model and Analytic Results 

83 We now investigate which instabilities drive the two distinct types of holes by examining 

84 in more detail the development of streaming instabilities. Using two drifting Maxwellians 
as to model the electron distribution and a single Maxwellian to model the ion distribution, 
ae we fit the distribution functions obtained from the simulations and substitute the theoret- 

87 ical fittings into the local dispersion function derived from kinetic theory for waves with 

88 < u < fi e [Che, 2009]: 

1 + + QZ(Q)] + 2(1 "^ [1 + /o(A)e- A Cei^(Cei)] (1) 

K u ti ft u tel 

+^[l + / (A)e- A C e2 Z(C e2 )] = 0, 

m where Q = (u - k z v di )/kv ti , Cei = (w - k z v del )/k z v ztel , ( e2 = (u - k z v de2 ) /k z v zte2 , 

90 A = k1v1 te /2Vl 2 e) 5 is the weight of the low velocity drifting Maxwellian, Z is the plasma 

91 dispersion function and I is the modified Bessel function of the first kind with order zero. 
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92 The thermal velocity of species j is defined by vfj = 2T t j/m,j and drift speed by i><#, which 

93 is parallel to the magnetic field (z direction). The electron temperature takes a different 

94 value along and across the magnetic field while the ions are taken to be isotropic. 

95 The fitting parameters of the distribution functions at Qtf = 3,4 are listed in Table 

96 1. The match between the parallel distribution and our fitted distribution is shown in 

97 Fig. 4 (a). We can see from the Table that the weight S of the low velocity electrons 

98 increases with time, indicating that momentum is transferred from the high velocity to 

99 the low velocity electrons. 

100 The theoretical 2D spectrum at f^t = 3 is shown in Fig. 4 (b). Two distinct modes 

101 are found, one with k parallel and the other with k nearly perpendicular to B. The peak 

102 of the parallel mode is around k z di ~ 20, which is close to the wavenumber of the cold 

103 plasma limit of the Buneman instability, k z di = u pe /vd e ~ 20. To confirm that the parallel 

104 mode is the Buneman instability rather than the two-stream instability, we exclude ions 
Ms from our calculations. The mode obtained only with electrons is shown in Fig. 4 (c). 
ioo The two-stream instability has a much smaller growth rate. Thus, the parallel mode 

107 is the Buneman instability. The peak of the nearly-perpendicular mode is centered at 

108 (k x di, k z di) = (22, 5). The frequency of this mode is ~ 13Qj which is in the LH frequency 

109 range for the present simulation so the nearly-perpendicular mode is the LHI [McMillan 
no and Cairns, 2006; Che et al, 2009]. 

in As a test of this interpretation, we compare the phase speed of the modeled waves across 

n2 (v px ) and along (v pz ) B with the simulation data. The assumption here is that since the 

us fraction of trapped electrons in any given electron hole is small, the non-trapped particles 
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n4 control the phase speeds of the wave and the linear dispersion characteristics can be used 

us to interpret hole propagation. It is well known that the Buneman instability can form 

U6 parallel bipolar structures. This instability, which has a very low parallel phase speed to 

n7 enable coupling to the ions, is the source of the electron holes moving slowly parallel to the 

us magnetic field. Thus, the LHI should be responsible for the oblique, fast-moving electron 

no holes marked by the downward arrows in E z and the oblique stripes in E x in Fig. 1. 

120 This interpretation is consistent with the parallel phase speeds v pz of the Buneman and 

121 LH instabilities obtained by the kinetic model which are shown in Fig. 2 (c). The phase 

122 speed of the Buneman instability with 9 ~ is close to zero. The three arrows from left 

123 to right (black,red and green) indicate the position 9 of the maximum-growing mode of 

124 the LH instability at f^t = 1, 3, 4 shown in Fig. 2. The phase speed of the LH instability 

125 is initially low and then increases to 4ca at Vtit = 3 and to 7ca at Vtit = 4 . The high 

126 phase speed of the LHI is consistent with the fast-moving electron holes seen at late time 

127 in the simulation. As a further check on this interpretation, in Fig. 5 (a) we stack the cuts 

128 of E x (x) along x at different times. The slope of the curves is the phase speed v px . We 

129 see that at Vtit = 3 v px ~ 0.6ca- In (b) is the theoretical phase speed v px at Qrf = 1,3,4 
wo calculated from the model. At fltf = 3 the v px of the LH wave, marked with the "*" , is 

131 around 0.6ca, consistent with the value from the simulation. 

4. Conclusion 

132 In summary, we have demonstrated through simulations and an analytic model that two 

133 distinct classes of electron holes are generated simultaneously in the intense current layers 

134 that form during magnetic reconnection. The sources of the holes are the Buneman and 



DRAFT 



January 19, 2010, 6:30am 



DRAFT 



CHE ET AL.: LH INSTABILITY AND ELECTRON HOLES X - 9 

135 LHI. The LH waves produce a transverse field E x as well as the bipolar structures E z that 

we trap electrons to form electron holes. These electron holes move along the magnetic field 

is? at the phase speed of the LH wave. Electron holes formed by the Buneman instability 

us move more slowly. The simultaneous existence of electron holes with two distinct phase 

us speeds enables electron scattering over a much larger range of velocity space than would 

140 be possible by either either instability alone. Electron dissipation in the intense current 

141 layers that form during reconnection is therefore enhanced. The LH electron hole was 
«2 also independently observed by 2D Vlasov simulations [Newman and Goldman, 2008]. 
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Figure 1. (a): The current sheet j ez in the x — y plane at fijt = 3.3. (b, c): The 
spatial structures of the electric fields E x and E z in the x — z plane in a cut through the 
current layer. 
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Figure 2. (a, b): Cuts of E x (z) and E z (z) around the x-line at different times from 
the simulation, (c): The theoretical parallel phase speed v pz vs. the angle 9 between 
wavevector k and magnetic field at flit = 1,3,4 (black solid, red dashed and green dash- 
dotted lines). The arrows denote the angle 9 of the fastest-growing mode of the LH 
instability at the three times in (d). (d) The theoretical growth rate 7 max of fast-growing 
mode vs. the angle 9 at the three times in (c). 
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Figure 3. (a): Spatial structure of E z at flit = 3 in the current layer, (b): The phase 



space z — v ez at x ~ 1.2 of (a) 
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Figure 4. (a): Electron and ion distribution functions f(v z ) around the x-line at 
Qit = 3 from simulations (blue solid) and the model (red-dashed) with the ion distribution 
function reduced by a factor of four. In (b) the 2D spectrum includes both electrons and 
ions and in (c) is without the ions. 
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Figure 5. (a): Cuts of E x (x) at different times from the simulation, (b): Theoretical 
phase speed v px vs. 6 at Q{t = 1,3,4, denoted by black solid, red-dashed and green 
dash-dotted lines. 
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